Report for Model 5

— for both jointly and separately

2024-10-08

1. Joint Modeling

1.1 Trace plots for convergence check

The current MCMC setting is:

  • 300000 iteration;
  • 200000 burn-in;
  • 100 thinning.

1.2 Gelman and Rubin’s convergence check

## Potential scale reduction factors:
## 
##         Point est. Upper C.I.
## HDevsum       1.01       1.04
## LDevsum       1.01       1.03
## dh0           2.26       4.47
## dh1           2.03       3.83
## dh2           1.15       1.55
## dl0           1.60       2.78
## dl1           1.44       2.53
## dl2           1.06       1.27
## 
## Multivariate psrf
## 
## 1.74

1.3 ACF Plots

Here we plotted ACF plots for the following variables:

  • Total deviance;
  • Variables that didn’t pass the convergence check.

1.4 WAIC results

LevelH LevelL
DIC 1223.23311 22091.4702
DIC3 1159.15231 22157.2817
PWAIC 43.55523 237.6256
WAIC 1186.24757 22178.4975

2. Separate Modeling of High-Level

2.1 Trace plots for convergence check

The current MCMC setting is:

  • 300000 iteration;
  • 200000 burn-in;
  • 100 thinning.

2.2 Gelman and Rubin’s convergence check

## Potential scale reduction factors:
## 
##         Point est. Upper C.I.
## HDevsum       1.03       1.13
## LDevsum       1.01       1.02
## dh0           2.51       5.01
## dh1           7.57      20.60
## dh2           4.89      15.02
## dl0           1.20       1.69
## dl1           1.16       1.55
## dl2           1.09       1.32
## 
## Multivariate psrf
## 
## 4.87

2.3 ACF Plots

Here we plotted ACF plots for the following variables:

  • Total deviance;
  • Variables that didn’t pass the convergence check.

2.4 WAIC results

H2
DIC 1375.85669
DIC3 1245.71574
PWAIC 90.38495
WAIC 1314.91290

3. Separate Modeling for Low-level

3.1 Trace plots for convergence check

The current MCMC setting is:

  • 300000 iteration;
  • 200000 burn-in;
  • 100 thinning.

3.2 Gelman and Rubin’s convergence check

## Potential scale reduction factors:
## 
##         Point est. Upper C.I.
## Devsum       0.999       1.00
## LDevsum      1.000       1.00
## dl0          1.001       1.01
## dl1          1.003       1.02
## dl2          1.003       1.02
## 
## Multivariate psrf
## 
## 1

3.3 ACF Plots

Here we plotted ACF plots for the following variables:

  • Total deviance;
  • Variables that didn’t pass the convergence check.

3.4 WAIC results

L2
DIC 22218.8518
DIC3 22261.9146
PWAIC 294.4887
WAIC 22293.6689